The induced charge in a Frohlich polaron: Sum rule and spatial extent. 



Sergio Ciuchi 1 Jose Lorenzana 2 'fcl and Carlo Pierleoni 1,2 

1 INFM and Dipartimento di Fisica, Universitd de L'Aquila, 

via Vetoio, 1-67100 Coppito-L'Aquila, Italy. 

2 INFM, UDR Roma I, Universitd di Roma "La Sapienza" , 

Piazzale A. Moro 2, 1-00185 Roma, Italy. 
(February 1, 2008) 

Within the path integral formalism, we derive exact expressions for correlation functions measuring 
the lattice charge induced by an electron and associated polarization in Frohlich polaron problem. 
We prove that a sum rule for the total induced charge, already obtained within approximated 
approaches is indeed exact. As a consequence the total induced charge is shown rigorously to be 
temperature independent. In addition we perform path integral Monte Carlo calculations of the 
correlation functions and we compare with variational results based on Feynman method. As the 
temperature increases the polaron radius decreases. On the other hand at high temperatures the 
electron motion is not hindered by the lattice. These apparently contradictory results are discussed. 
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I. INTRODUCTION 

An electron added to an insulating polar crystal foci 
a quasiparticle called dielectric polaron after FrohlichM 
This has been recognized as a fundamental field theoret- 
ical problemEI. More recently a variety of novel mate- 
rials have emerged which present interesting properties 
when doped away from an insulating phase, like colos- 
sal magnetoresistent manganites and high temperature 
superconducting cuprates. The fact that these are po- 
lar crystals has produced a renewed interest in Frohlich 
polaron problemm! 

Roughly speaking a dielectric polaron is composed of 
an electron and the (opposite) charge that it induces in 
the lattice. Electron and induced charge attract each 
other so that for the electron to move it has to drag the 
induced charge resulting in an increase of the quasipar- 
ticle mass. In this work we study correlation functions 
which measure the magnitude and spatial extent of the 
induced charge and associated polarization field. 

We derive a rigorous sum rule which states that the 
total induced charge equals the charge induced by a clas- 
sical (static) electron and it is independent of tempera- 
ture. This result, welLknown within perturbativecl and 
variational approaches^, is proven here to be exact. The 
large distance behavior of the electric field is determined 
by this sum rule. In addition we discuss the short dis- 
tance and the high temperature asymptotic limit of these 
quantities. These results provide constrains to approxi- 
mations on the polaron problem. 

Real space path integral Monte Carlo (PIMC) method 
is used to evaluate correlation functions. These are 
compared, with Feynman's variational approximation 
(FVA)li!ffl and analytical results in weak and strong cou- 
pling. We find that the polaron radius is determined 
at low temperatures by the electron-phonon coupling a 
alone while at high temperatures it is proportional to 



the de Broglie thermal wave length (At = h^2irP/m 
with (3 the inverse temperature) and becomes indepen- 
dent of coupling. We define a polaron crossover temper- 
ature T*(a). Although the electron induces a tempera- 
ture independent charge in the lattice the induced charge 
hinders the electron motion only below T*(a). At high 
temperatures thermal effects wash out the hindering ef- 
fect of the induced charge but, we remark again, not the 
induced charge itself. 



II. ANALYTICAL RESULTS 

Our starting point is the effective action for the 
Frohlich polaron problem, after the phonons have, been 
eliminated with the usual path integral techniquesu, S = 
Sq + Si with 
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here m is the electron mass, a — e 2 m}/ 2 /TiI{2Tiujl) 1 ^ 2 , 
is the coupling constant with uil the phonon frequency, 
1/e = I/cqo — 1/eo and e w (eo) is the high (zero) fre- 
quency dielectric constant. 
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is the phonon propagator. 

We are interested in the correlation function between 
the electron charge density — en(r) and the charge in- 
duced in the lattice eni(r) normalized to the probability 
density to find an electron at a given point i.e. the inverse 
of the volume V. Dropping the charges this is defined as, 
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(n(0)ni(r))/y- x = g(r)/e 
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Averages are defined as path integrals weighted by S 
_ /Px e - s / f '(...) 
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where the paths entering in Eq. (Q) depart and arrive at 
the same point. Further integration over such a point is 
not performed and this assigns to Eq. (^) the meaning of 
a constrained average with x(0) ~ x(/3) = 0. Those av- 
erages are however equivalent to the unconstrained ones 
because of translational invariance. We used the spher- 
ical symmetry of the problem to define g(r), and we di- 
vided by I on the right hand side of Eq. (||) for later 
convenience. Another quantity of interest is the induced 
lattice polarization 



P(r) = (n(0)P(r))/V- 



(5) 



VP(r) = — eni(r) is the density of polarization opera- 
tor. The correlation function in Eq. (5|) is related to the 
induced electrostat ic p otential [W 2 V(r) — —4:Treg(r)/e] 
considered in Ref. [Il] and which will not be discussed 
here. We stress that these quantities have the meaning of 
correlation functions measuring average induced charge, 
polarization and induced potential at a distance r from 
the electron position. _ 
The charge density operator for the phonons isE3, 



kLJL ]>>Q k e ikr . 
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where Qk is the dimensionless displacement for momen- 
tum k phonons. Inserting Eq. (||) in Eq. (||) we obtain 
an equation for g(r) as a function of the density dis- 
placement correlation function (n(O)Qk)- The phonon 
variables can be traced out by standard methods. We 
have found that it is possible to give an exact expression 
for the correlation functions in terms of path integrals 
weighted by the effective electronic action of Eq. (|l|) . We 
find for the density-induced-density correlation function 



g(r)= / drU(T)(S[r-^T)}) 
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and for the polarization field 



P(r) 
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where f = r jr and 



sinh(cj L T) + sinh[cj L (n/3 - r)\ 
1 ' L 2tanh(/3?iw L /2)sinh(/3?iw L ) ' 1 ' 

Within FVA, the variational quadratic action can be ex- 
ploited in Eqs. to analytically perform the aver- 
ages and to recover the results of Refs. [lO] and O. 



Eqs. (0),(||) have a simple physical interpretation. The 
induced charge can be seen as "distributed" along the 
electron path with weight U (r) . The polarization is the 
superposition of polarizations associated with those ele- 
mentary induced charges. 

Eq. (0) can be integrated in the whole space using the 
properties of the Dirac's S function. Since / (1tU{t) = 1, 
we conclude that g(r) is normalized to one. The total 
induced charge q is computed by integrating the density- 
induced density correlation function in Eq. (H): 
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which completes the proof of the sum rule. The total in- 
duced charge amounts to the charge the electron would 
induce if it where a static classical particle. In other 
words there are no quantum corrections to the total in- 
duced charge. 

An alternative derivation can be work out following 

QuemeraisEl. From the time derivative iKP = \P,H] 
one obtains 



(11) 



Here H is the Hamiltonian, D is the electric displacement 
operator due to the electron (V.D = —Aiten) and iKP = 
[P, H]. Taking the divergence we obtain a relation for 
the charge operators: 



V.P(r) 
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We can integrate this expression in all space and take the 
thermodynamic average. The left hand side is propor- 
tional to the average of the net force felt by the lattice 
at the boundary of the system which should vanish at 
equilibrium. The right hand side gives Eq. ([To|). 

Eq. (|l0|) shows that the induced charge is independent 
of temperature. This contradicts the naive argument that 
all polaron effects should disappear at high temperatures. 
To understand this behavior one can do an analogy with 
the behavior of an harmonic oscillator in an external field. 
In that case, because of harmonicity, one gets a displace- 
ment which is temperature independent. Here roughly 
speaking the harmonic oscillator represents the phonon 
coordinates and the external field is the field produced 
by the electron on the phonons. The induced charge is 
a measure of how much the ions displace from their bare 
equilibrium positions in the presence of the electron. As 
for the single harmonic oscillator, this "displacement" is 
independent of temperature. Only anharmonicities can 
make the induced charge temperature dependent. 

Using the sum rule it is easy to see that at distances 
much larger than the polaron radius, as defined be- 
low, the polarization field goes as P(r) = — er/(Awer 2 ). 
Clearly the distortion produced by the electron is long 
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range, a^iact that is not always recognized in the 
literaturecJ. The total electric field (always in the sense 
of a correlation function) is given by E = D — AirP where 
we should include in D = — ef/(e oc r 2 ) the high frequency 
screening. At long distances we have E(r) = — ef/(eor 2 ) 
which means that the electric field generated by the elec- 
tron gets screened by the static dielectric constant. This 
is generally expected but to the best of our knowledge 
has never been proven for all coupling and temperatures. 

Now we discuss the short distance behavior. At dis- 
tances much smaller than the polaron radius we expect 
that the effect of the interaction becomes irrelevant in 
the functional integrals. This is because the latter are 
dominated by electron paths with short wave length or 
equivalently high kinetic energy. We can then replace the 
total action by the free electron action in Eqs. We 
obtain the asymptotic result: 



lim P(r) 

r— vO 
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where I = yJh/2mujL is the harmonic oscillator charac- 
teristic length. Using the same argument we obtain that 
g(r) oc r _1 for r — > and the proportionality coeffi- 
cient can also be obtained with the same method. The 
latter behavior has been obtained in Ref. ^ within the 
FVA. These results coincide with lowest order perturba- 
tion theory. 

At high temperatures we also expect that the effect of 
the interaction becomes irrelevant and so we can replace 
again the total action by the free electron action in the 
functional integrals. The high temperature asymptotic 
result for g(r) is 

^ 9{r) - vqki exp (p^r) phuJL 1 (14) 

This result has also been obtained in Ref. ^ within FVA. 
We remark that although the density-induced density 
correlation function does not vanish for large tempera- 
tures the polaron effective mass tends to the bare electron 
mass. 



III. NUMERICAL RESULTS 

Now we discuss the spatial extent of the induced charge 
at general couplings and temperatures. We have evalu- 
ated averages in Eq. |]) using PIMC. Eqs. (0),(§) being 
expressed in real space rather that in Fourier components 
are more suitable for this purpose. We have performed 
Metropolis PIMC calculations within the imaginary time 
discretization scheme. In order to regularize the attrac- 
tive divergence of the retarded action at short distance, 
and to improve the convergence with the number of the 
imaginary time discretization points, we have developed 
a preavpeaging procedure similar to the one used for local 
actionsE£lc3. Details of the method and more extensive 



results will be given in a separate publication^ Here we 
just say that the results for the g{r) are well converged as 
checked by doubling the number of imaginary time slices. 
Previous MC studies of the Frohlich polaron -were limited 
to small (a < 4)c3 or intermediate (a < 7)113 couplings 
and wece-ipcused to the calculation of the ground state 
energyMLa and effective masst£l. 

Following Ref. ^ we have also computed g(r) in the 
FVA, i.e. using Feynman's quadratic action to evaluate 
the average appearing in Eq. (0)E£|. 

In Fig. |l| we show 47rr 2 <?(r) for weak, intermediate and 
strong coupling. We also show the weak coupling re- 
sult obtained by perturbation theory at T — and the 
strong coupling result in the Landau-Pekar approxima- 
tion. For all couplings the correlation function decays 
exponentially with distance as expected for a polaron. 
The area under the curves is one according to the sum 
rule Eq. (|l0|). The short distance asymptotic behavior 
Eq. (|RJ) is exactly satisfied in FVA and is also satisfied 
in the PIMC within the numerical error. We define the 
polaron radius r m as the distance at which Anr 2 g(r) is 
maximum. As the coupling increases the polaron shrinks 
indicating the progressively more localized nature. 




FIG. 1. 4nr g(r) as a function of r for different couplings 
and approximations. From left to right bold curves (FVA 
after the formalism of Ref. ^) and points (PIMC) corresponds 
to a — 12,6,1, respectively, and f3hu)L ~ 20. Thin solid 
line is perturbation theory, thin dashed line the Landau-Pekar 
strong coupling approximation both at PTiujl = oo. 

In Fig. U we show the temperature dependence of 
4irr 2 g(r) at intermediate coupling. 

In Figs. UU we see that, apart from a small overesti- 
mate of the polaron radius, FVA reproduces fairly well 
the PIMC data. We can then safely use the FVA to dis- 
cuss the spatial properties of the polaron. Notice that 
the good agreement found between PIMC and FVA is 
not obvious since in principle only the free energy is ex- 
pected to be accurate for the latter. 

Contrary to the naive expectation the effect of tem- 
perature is to shrink the polaron (Fig. ^|). Our phys- 
ical explanation is the following: At low temperatures 
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phonons relax shrinking the spatial extent of the elec- 
tron till the increase in electron kinetic energy balances 
the gain in electron-phonon interaction energy. At high 
temperatures, however, a typical electron has energy 
E = 3/2/3 and momentum hk = \j7>ml (3. One can con- 
struct a wave packet of width Ak in momentum space 
using plane waves with higher and smaller energy with- 
out affecting the electron internal energy. The biggest 
Afc which will not affect the electron internal energy is 
of order of k itself. One can then achieve a localization 
of the electron of order 1/fc ~ h/^3m/(3 ~ At- It fol- 
lows that the phonons can relax at practically no cost 
till the polaron radius stabilizes at a value of this or- 
der. In fact at high temperatures the asymptotic value 
of the polaron radius can be obtained from Eq. (14): 
lim^o r m = W '(3fruJL /p-i= 0.2Ar- This scaling has been 
found by Sethia et alE3 for the mean square displace- 
ment of the electron in imaginary time. Notice that the 
polaron radius becomes independent of the coupling. 

The same behavior of r m has been obtained in Ref. |ll| 
within FVA. The authors of Ref. |ll| ascribe the high tem- 
perature behavior of the polaron radius to the increased 
fluctuations of the phonon field. We conversely think that 
the thermal fluctuations of the electron are responsible 
of the high temperature behavior of the polaron radius 
and the phonon field acts only as a probe of the intrinsic 
thermal length of the electron as explained above. 




FIG. 2. 4irr 2 g(r) as a function of r at a = 6. From left to 
right curves (FVA after the formalism of Ref. ^ ) and points 
(PIMC) refers to (5Tiuj l = 0.1, 1.0,20, respectively. 

To characterize the temperature and coupling depen- 
dence of the polaron size we plot in Fig. || the polaron 
radius as a function of coupling for different inverse tem- 
peratures. We see that when the temperature is such 
that 0.2At > I (low temperatures, see the curves for 
(3hu>L = 20, oo), the polaron radius exhibits little temper- 
ature dependence at all couplings and simply interpolates 
between the weak coupling polaron radius r m — I and the 
Landau-Pekar polaron radius r rn = 3/^/ n/2/a. When 
0.2At < I two different regimes occur. At small coupling 
the polaron radius tends to saturate at 0.2At (horizontal 



arrow) whereas at high couplings one recovers the low 
temperature polaron radius r m (a,T = 0). We can define 
a crossover line when these two lengths are of equal mag- 
nitude so that the crossover temperature as a function of 
a is given by the equation r m (a, T = 0) = 0.2At(T*). In 
the FVA approximation^ we find the crossover temper- 
ature in energy units to be T*(a) — Q.15hu}Lv(a,T = 0) 
with v the usual Feynman variational parameter!! For 
a — > 0, T* goes to a constant of order of O.BTllul whereas 
for large a it increases quadratically with a. 
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FIG. 3. The polaron radius as a function of coupling for dif- 
ferent inverse temperatures. Results from FVA are shown as 
dash-dotted (PTilol = 0.1), dashed (J3Jujl = 1.0), thick-solid 
(/3huJL = 20) and thin-solid lines (JSTiwl = oo). Dotted line 
is the Landau Pekar approximation. Results from PIMC 
calculations are shown as triangles (J3Tuol = 0.1), circles 
(f3hu)L ~ 1.0) and squares (PTujjl = 20). The horizontal ar- 
rows indicate the value 0.2\t/1 for f3Tiu)L = 1 (upper) and 
I3Tiu)l =0.1 (lower). The arrows for /3hu>L ~ 20, oo are out of 
scale. 

In the low temperature regime (T < T*) the polaron 
radius becomes almost temperature independent and is 
determined by the couplingE3. The high temperature 
regime T > T* is characterized by a polaron radius which 
is independent of the coupling and is determined by the 
temperature alone (r rn — 0.2At)- 



IV. CONCLUSIONS 

We have studied the charge induced in the lattice by an 
electron in Frohlich polaron problem and the associated 
polarization field. We have derived relations which ex- 
press the charge-induced density correlation function in 
terms of path integral involving only the electronic degree 
of freedom which are suitable to be evaluated by PIMC 
method. A rigorous sum rule was derived that determines 
the total induced charge and the long distance behavior 
of the polarization field. We give also the asymptotic 
limits of these quantities at short distances and at high 
temperatures. We have compared results obtained using 
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FVA through the lines of Refs. |9|-[Tl| and those obtained 
by PIMC method. To the best of our knowledge this 
is the first PIMC computation of real space correlation 
functions in Frohlich model. From the spatial depen- 
dence of the induced charge we obtained a polaron ra- 
dius. The polaron radius is determined by the coupling 
at low temperatures and by the thermal wave length at 
high temperatures with a crossover temperature that we 
evaluated in FVA. 

At high temperature a polaron with small radius and 
small effective mass is achieved. These results are not in 
contradiction because the small radius at high temper- 
atures is a thermal effect of the electron and it is not 
related to the lattice response. The lattice acts only 
as a probe of the intrinsic electron localization radius 
namely At- Obviously this small radius polaron has 
nothing to do with the Holstein zero-temperature small 
polaron which induces almost local lattice displacements 
and moves coherently with a large effective mass. 

The PIMC polaron radius is always smaller than the 
FVA calculation in the range of coupling and tempera- 
ture studied. This effect is more pronounced at inter- 
mediate couplings. The overall temperature dependence 
agrees with the findings of Ref . p|-pd| however our physical 
interpretation is different. 
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